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Abstract 

The DO collaboration has announced the observation of the like-sign dimuon charge asymmetry 
since 2010, which has more than 3a deviation from the Standard Model prediction. One of the 
promising explanation is considering the existence of flavor changing Z' couplings to the b and 
s quarks which can contribute to the off-diagonal decay width in the B s — B s mixing. Model 
construction is highly constrained by the recent LHCb data of lfb -1 integrated luminosity . In 
this paper, we analyze the experimental constraints in constructing new physics models to explain 
the dimuon charge asymmetry from the CP violation of the B s system. We present limits on Z' 
couplings and show that it is impossible to obtain the la range of the dimuon charge asymmetry 
without the new contribution in the B^ system. Even with arbitrary contribution in the B^ system, 
the new couplings must be in the fine tuned region. 
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I. INTRODUCTION 



The like-sign dimuon charge asymmetry from the semi-leptonic (s£) decay of B s ^ meson 
is given by, 

A b - N++ ~ N ~~ (1) 
ai ~ N++ + N~ ' { ) 

where N ++ corresponds to each B hadron decaying semi-leptonically to fi + X, and similarly 
N to fi~X. In 2010, the DO collaboration at the Tevatron announced the first observation 
of the large dimuon charge asymmetry, which deviated about 3.2a from what is expected 
in the Standard Model (SM) [lj. In 2011, the result from the analysis with 9 fb _1 data was 



announced as 

A b s e = -(7.87 ± 1.72 ± 0.93) x 10~ 3 , (2) 
which has about 3.9cr deviation from the SM prediction 

A b * M = (-2.8i°;jj) x 10- 4 . (3) 

To explain the observed asymmetry, we need additional sources of CP violation from the 
new physics (NP) beyond the SM in the B Sjd mixing and/or decay. 

The contribution from each neutral B° and B° s meson is parametrized by the flavor specific 
asymmetry 

ad _ T(B d ^fi + X)-T(B d ^fi-X) 
a ° l ~ T(B d -> fi+X) + T(B d -> fi-X) ' 
T(B S -> fjL+X) - T(B S -)• ii- X) 



a% = 



T(B S -> fi+X) + T(B S -)• n~X) 



The fraction of each flavor specific asymmetry in the total asymmetry A b si at the Tevatron 
energy 1.96 TeV depends on the mean mixing probabilities and the production fractions of 
B° and B° s mesons such that 

A b si = (0.594 ± 0.022)a^ + (0.406 ± 0.022)t& , (5) 
which leads to 6 : 4 production of the like-sign dimuons from the bd(db) and bs(sb) mesonsJl] 



1 This is different from the 2010 prediction of about 5 : 5 production. 
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Imposing the lower limits of the muon impact parameter (IP), it is possible to reduce the 
background dramatically, which is mainly from the long-lived charged mother particles of 
the muon and the anti-muon. In the 2011 data, the separation of the sample by the muon 
impact parameter provides the separate determination of a d e and a s sl such that 

a 1 ^ = -(18.1 ± 10.6) x 1(T 3 , (6) 
a d si = -(1.2 ± 5.2) x 10" 3 , (7) 



where the SM predictions using the SM fit of \V ub \ = (3.56+°;2o) X 10" 3 ^ are 

a^ SM = (1.9 ± 0.3) x 10" 5 , (8) 

4 SM = -(4.1 ±0.6) x HT 4 . (9) 

The separately determined a s sl has about 1.7a deviation from the SM prediction if a d e can 
be freely chosen to fit the data. Similarly the a d st is within la if a 8 st can be arbitrary. It 
should be noted however that in order A b s£ to be within la from its measured value, a large 
contribution from new physics in a s s£ is necessary as we see bellow. 

We reproduced the x 2_ fit combining the impact parameter cut (120/im) analysis of IP<i2o 
and IP>i2o m Fig- CD As seen in the figure, we need to consider both of a d sl and a s si together 
in the 2D plane. In this figure, we used the central values in the fraction of contribution 
by a d e and a s sl in A b s£ , shown in In Fig. CD (b), the red dot denotes the observed 
central values of (a d e , a s si ) while the black dot represents the SM predictions. We chose three 
representative values for a s si and one for a d e , whose distances form the SM values are given 
in terms necessary enhancement of — anc ^ a f^/( a i) SM - First of all, if we allow 

arbitrarily new physics in a d e , we still need at least 68 times bigger size of a s s£ to explain the 
asymmetry within la where g^/(g^) sm = 21 at the la boundary. Without NP contribution 
to a d £ , the enhancement — a s s£ / (a s s£ ) SM > 580 is needed to explain the asymmetry within la. 
To be at the central point in which the \ 2 fit is the best, a s sl should be enhanced by a factor 
of 950 while small enhancement in a d e /(a d e ) SM < 3 is enough. Therefore, we conclude that 
the observed value of A b s£ requires large NP contribution in a^O We take three points as 
references. 

(4/(^) SM , ayM SM ) = (21, -68), (1, -580), (1, -950) , (10) 



2 If we aim for the asymmetry within 90% confidence region (1.65cr) or 2a region, the observation result 
can be achieved only by the enhancement in ao without having any contribution in a s st . 
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FIG. 1. We reproduced the measurements with different muon impact parameter (IP) selections 
according to [2] in (a). The bands are the 90% uncertainties on each individual measurement 
of IP<i20 (Gray), IP>i20 (Cyan), and the result without the IP cut (Purple) in Eq.([2]). The 
green (68%), yellow (90%), and orange (95%) ellipses are obtained from the x 2_ fit combining the 
measurements of IP<i20 an d IP>i20 using the independent data sample. In FigQ] (b), the red dot 
denotes the observed central values of (a^ e ,a s se ) while a black dot represents the SM predictions. 
We chose three representative values for a s sl and one for eA, whose distances form the SM values 
are given in terms necessary enhancement of — o, s se / (a s s£ ) and o^e/ ( a si) • 

where the small enhancement in a d st at the third reference point is not considered. 

As a promising example explaining the large dimuon charge asymmetry, the Z' scenario 
with both flavor diagonal and off-diagonal couplings has been analyzed j^, 5|. In this paper, 
we study the validity of Z' boson explanation satisfying Eq. ffTOj) by checking recently updated 
experimental constraints from the B/B s meson decays and mixing Especially, the recent 
LHCb results provide very strong bounds. We study the operators (sx7 m &x')('7V7m 7 v) or 
(sx7 M &x)(cy7 M cy) (where X, Y — L,R), because NP contribution to b — > st + t~ is weakly 
constrained from the Bi(B s — > r + r _ ) [8J, and the effect on a s sl from (sx7 M &x)(cy7 At cy) can 
be enhanced by the interference with the W boson exchange [5|. Unlike other papers, we 
present our results in terms of the actual Z' couplings for a fixed Z' mass, Mz>- (Readers 
3 The electroweakprecision test results can also provide strong constraints when the mixing of the Z' and 



Z boson exists 



7[, while we do not consider such effect here. 
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can simply rescale constraints on the couplings for different value of Mz>-) Therefore, it will 
be easy to see the feasibility of realizing allowed space of Z' couplings from the view point 
of model building that we don't discuss in this paper. The effective set-up only considering 
the Z'sb and Z't + t~ [Z'cc) couplings is used, regardless of their theoretical origins. 

While we were in the completion of our work, a similar analysis for the operator of 
(sx7 M ^x)(cy7 A tCy) was appeared [9|. They chose special cases either one of the couplings of 
Z'clCl and Z'crCr is turned off or they are set to equal. Comparing to this simplification, 
our analysis deals with general case with more systematic approach. By doing this, we point 
out the Z'cc couplings must be (almost) axial vector-like from the constraint B s — )■ J ftp <fr 
and quantitatively see how much the axial vector relation can be violated by combining 
other experimental bounds. For the Z't + t~ couplings, we note that our analysis includes 
the constraint from b — > svv that has not been discussed in the preceding studies. 

This paper is organized as follows. We provide a summary review of the Z' explanation 
on the dimuon charge asymmetry in Sec. [Til Then, we analyze the current experimental 
bounds constraining the NP model construction explaining the asymmetry and apply the 
bounds to the Z' properties in Sec. IHII The experimental results we will analyze contain the 
measurements of the mass difference AM S and the width difference Ar s after the mixing. 
We also included the bounds from the CP violating phase (f) J J^^ of the B s — )■ J ' j 'ip process, 
the inclusive b — > svv, and the sin 2/3 from the golden plate B J/ipK s . In Sec. HVlandlVl 
we directly obtain the combined constraint on the Z' model parameters in the models with 
the Z't + t~ coupling and the Z'cc coupling, respectively. Finally, we give the conclusions in 
Sec. ED 



II. THE LIKE-SIGN DIMUON CHARGE ASYMMETRY 

The B q — B q oscillations for q = s,d are described by a Schrodinger equation 

(11) 





where M and T are the 2x2 Hermitian mass and decay matrices, which are dispersive and 
absorptive parts in the time dependent mixing respectively. The differences of masses and 
widths of the physical eigenstates are given bv the o ff -d lag o„a, events a8 Q 

AM, = 2|M 1 9 2 | , Ar g = 2|ry cos(p q , (12) 



up to numerically irrelevant corrections of order m^/M^ as long as AM 3> Ar for B q meson 
system. The CP phase difference between these quantities is defined as 

g = Arg. , (13) 



r 

1 12 



where the SM contribution to this angle is 



0SM = ( _ 7>5 ± 2 A) x 1Q -2 ^ ^SM = (3>8 ± l V) x 10 -3 _ (14) 

The flavor specific charge asymmetry a q st is related to the mass and width differences in 
the B q — B q system as 

\r q I Ar 
• ; A*"' "35, 



•« = Im iif = iiit| Sin0 « = AJit tan ^' (15) 



iere, the experimental value of AM S obtained from the LHCb 0.34fb 1 with 68.3% CL. is 



AM S = 17.725 ± 0.041(stat.) ± 0.026(sys.) ps" 1 , (16) 
(the combined result of CDF and DO is AM, = 17.78 ± 0.12 ps _1 ) while the SM prediction 



is 



12| 



(AM S ) SM = (17.3 ±2.6) ps" 1 (17) 
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which corresponds to f Bs =231 MeV and B B = 1.28 of Eqs. j22J and ()23]) |3|. 

The observed value of AM S has not so much deviated from the SM prediction. Without 
considering the NP contribution to Tl 2 , therefore, it is impossible to obtain the observed 
central value of a s st from Eqs. ( fl5|) and ( TlBl) for q = s even we assume sm<f) s = — 1. With 
the recent LHCb bound for <$> J J^^ , the maximally possible enhancement of a s sl in this case 
is outside the boundary of la of the observed value in (flOj) as seen in Fig. |3] (b). Therefore, 
an additional NP contribution to T\ 2 is preferred to explain the like-sign dimuon charge 
asymmetry through the B s — B s mixing. 

To probe the NP contribution, we split Ti2 or M\2 to the SM and NP contributions as 



for real and non-negative parameters h q and h g , with the phases constrained in the region, 
< a q , a q < 7r. Then, the flavor specific charge asymmetry is given by J4] 

a I 1 12 l 



|M 1 9 9 SM | 1 + h% + 2h q cos 2a q 



-h sm2a (l + h a cos2cr ) + h„ sin2<j (l + h„ cos2<5" ) \ cos 



lSM 



+ j(l + /i 9 cos 2<Tq)(l + h q cos2cr g ) + sin 2a q sin 2cr g | sin ( 



N SM 



(19) 



Also, the ratio of Eq. (jT9l) to its SM value is given by 

1 



-a q J(d 



1 + hf, + 2h q cos2cr g 
x i h q sai2a q (l + h q cos2a q ) — h q sin 2<j g (l + /iqCOs25"g)| cot0^ M 

(20) 

^1 + /i g cos25-q)(l + h q COS 2(T 9 ) + /i^/ig sin 2a q sin25"q| 

Note that the factor l/(l + h 2 q + 2h q cos2(x 9 ) in Eq. ((20]) is fixed by the ratio of AM^ M /AM q , 
near to 1. Therefore a sizable NP contribution to |r'^ p /r^| M | = h q is necessary if can take 
the dominant role in explaining the observed dimuon charge asymmetry. 

The Z' models to enhance the a s sl require the existence of nonzero off-diagonal couplings 
g^ b and gf b , where is the coupling of Z' to fermions ipL,R and Xl,r- Turning off one of 
the couplings g^ b and gf b for simplicity, this scenario demands the existence of rather large 
couplings \g^ R \ > 1 to explain the asymmetry within la range from the observed central 
value due to the strict AM S constraint. The situation is the same even in the case that 
the mass of Z' is similar to that of the Z boson. Such large g^ R couplings can violate the 
observations in the electroweak precision test (EWPT). Therefore, we need to turn on both 
of the flavor changing couplings g^ b and gf b . The scenario considering the g^ R couplings to 
explain the dimuon charge asymmetry will be called as Ll g TT scenario" in this paper. 

On the other hand, considering the nonzero Z' coupling to the charm quark pair can also 
explain the dimuon charge asymmetry by considering the interference of the NP contribution 
and the SM process. Due to the interference, the couplings g^ b R gcc R contribute to a s s£ linearly 
while g^' b R g^ R do quadratically so that the interference effect dominates the enhancement 
of a s si unless the NP contribution is larger than that of the SM. Therefore, it is possible to 
explain the asymmetry with rather smaller Z' couplings in this scenario so that we can avoid 
the direct constraint such as the decay of B s — > DD S j^]. The scenario considering such 



contribution will be called as u g cc scenario" in this paper. Describing the corresponding 
in each of our Z' scenario, there are six real free parameters, i.e. the complex g^ R and the 
real g^ R (<?,^' R ) since the diagonal couplings have to be real. 

Every experimental result depends not only on the mass of Z' but also on its couplings 
to the matter because the new interaction depends on the ratio (g^' R /gi)(Mz/Mz'), where 
gi = g/ cos 9w for g is the SU(2)x coupling and 6w is the weak mixing angle. Therefore, the 
experimental bounds can be applied for any values of Mz> by proper rescaling of the couplings 
gffc- According to this fact, we set the reference value Mz> = Mz for the representation 
of our analyses so that one can easily see the results for any Mz' one wants to analyze, by 
simple rescaling of the Z' couplings. Actually, our reference value of Mz> is not unrealistic 
since the 6-quark forward-backward asymmetry A b FB at the LEP can be explained in terms of 
Z' where Mz< ~ Mz and the non-zero g^ R and g^ R exist [3]. As a conservative approach, 
one can consider the heavy Z' whose mass is much larger than 1 TeV to avoid the current 
experimental limits when the Z' couplings to matter are SM-like By simply rescaling 

our final result in such a case, some Z' couplings to the matter should be much larger than 1 
to explain the dimuon charge asymmetry, which is unrealistic in the perturbative regime. On 
the other hand, one can also consider very light Z' cases whose couplings are small enough to 
avoid the direct Z' search bounds. Then, one needs to apply the other experimental bounds 
which we will explain from now on. 

The NP models accommodating the sizable new contribution in r\ 2 suffer from the various 
experimental bounds, mainly due to the recently updated LHCb data of lfb -1 . In the next 
section, we analyze the related bounds in detail by focussing on the enhancement of a s si . 

III. EXPERIMENTAL CONSTRAINTS 

In this section, we analyze the various experimental constraints in obtaining the new 
sizable contribution to r{ 2 from the B s — B s mixing. The NP contribution to rf 2 via the 
operator (sb)(ff) where / is a SM fermion can also affect the various B s or B meson 
decay processes □. In the Z' models, the new contribution is realized by the tree level 
FCNC process, which can be large enough to threaten the current experimental bounds. In 
this section, we introduce the experimental constraints from AM S , AT S} (j) J J^^, b — > suu, 



4 B generically denotes B®,Bf mesons. 



S 



and B — > J/ipKs- Then, we will show what extent the NP parameter space explaining the 
dimuon charge asymmetry can be constrained by such bounds, by applying our Z' scenarios. 

For the simplicity, we turn off the couplings g e l for the light leptons £ = e~ , /i~ not 
to consider the tree level NP contribution in the observations such as B — > X s £ + £~, B — > 

n l r 

K*£ + £~ } and B s — > £ + £~ as shown in |4|. For the case g bb ' ^ 0, a one-loop induced NP 
contribution can affect the b — > s^y. Such contribution is well summarized in our Appendix 
lAI for the future use. 



A. AM, 

The experimental measurements of AM S both from the LHCb and the Tevatron have 
no significant deviation from the SM prediction. Therefore, the allowed parameter space is 
highly constrained as shown in Fig. [2] in terms of the general parameters h s and 2a s . 

In the SM, the dominant contributions to M 12 come from the top quark loops and their 
effects are summarized as follow. 

Mf 2 M = ^ 2 Ml{V tb V:fS Q {ml/Ml)m B J 2 BsV2B B B (21) 



Here, Sq(x) is an Inami-Lim function for the corresponding box diagrams [16j, and t] 2 b and 
Bb are fj,b andujy independent quantities at a given order of QCD corrections. At the NLO, 
r]2B — 0.551 ll7| and Bb is given as 

B B = M/^)r 6/23 (l + BX LL ^ h ) , (22) 

where J5 = 1.627 (in NDR and / = 5) [17] . B\ LL is a bag parameter of a matrix element 

(B s \OX LL {n)\B s ) = \m\f B BX L \y) , (23) 

where 

0\ LL = {s L Yb L ){s Ll ^b L ) (24) 

and mg s and fg a are B s meson mass and its decay constant, respectively. 

For the evaluation of Eq.flm we use G F = 1.16637(1) x 1(T 5 GeV" 2 , M w = 80.399(23) 



GeV, m Bs = (5366.3 ±0.6) MeV, \V tb \ = 0.999152^™", \Vu\ = (4-03^) 18|. |] For 



5 Here and after, the figures in parentheses after the values give the 1-standard-deviation uncertainties in 
the last digits 
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FIG. 2. The yellow colored region denotes the parameter space allowed by the 90% C.L. (1.65cr) 
experimental bounds of AM S observed at the LHCb 0.34fb -1 . The rough upper limit of h s is 2.3 
according to this figure. The results from the CDF and DO are not so much different from this. 



the top-quark mass, we use m? t ole = 173.2 ± 0.9 GeV (correspondingly rnt(m f ) = 165.8 ± 0.9 
GeV ) [19j. Finally we use f Ba = (229 ± 6) MeV and B B = 1.291 ± 0.043 [20]. 

For the Z' and its flavor violating interactions of g^sij^biZ'^ and g^SR^bnZ'^, following 
effective operators, in addition to 0\ LL ', are induced at the scale where the Z' is integrated 
out. 



OX RR = {s R Yb R ){s Rl ,b R ) 

0{ R = {s L rbL){s R lM 
RL = (s R Yb R )(s Llll b L ) 



(25) 
(26) 
(27) 



At the same time, QCD corrections to the operators of Eqs. (1261) and (1271) induce following 
operators as well. 

O hR = (s R b L )(s L b R ) (28) 

RL = (s L b R )(s R b L ) (29) 

Using those notations, we write down an effective Hamiltonian at the scale fib as 
1 



n/Z' 
H eff 



2M% L 



V LL m 2 + (9 R b ) 2 ) OX LL + 2 <fa& + 2 rj^g* O l 2 r 



(30) 
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Note that we identify Of Y with Oj x (X, Y = L or R) at this stage, reflecting the fact that 
the QCD is vector-like and corresponding matrix elements are equal. The QCD corrections 
are given at the NLO in Ref. 2JJ as 

LL 6/23 . Oi s {m b ) ( 6 / 23 \ 

' + 4?r (^1-63 (1 — 775) ?7 5 / J 



rj 



3/23 . 
% + 



a s (m b ) 



An 



0.93% 24/23 + t? 3 / 23 (-2.10 + l.rfrfe)) 
((-11.73 + 0.78r7 5 )r7, 



3/23 - 24/23 \ . 0! s (m b ) 



3/23 + ^-24/23 3Q + ^25^ 



where ^5 = ai 5 \/J,z)/ci.s 0} (^b)- Parametrizing the hadronic matrix elements as 



.,(5)/ 



(B s \O l 2 r {^)\B s ) = 1 



3 \ m & + m i 



2 \m b + m s 

we get following expression for the Z' contrition to M 12 



(31) 
(32) 



M, 



12 



6M|, 



^ L ((^) 2 + (4) 2 )5i VLL -^ 



sb 



nib + ^ 



LR D LR 



'III D \ 



° LR R LR 
2 ''21 £> 2 



(33) 



For the evaluation of Eq.( )33|) . we use the two- loop RG running with the input of a s (Mz) = 
0.1184(7), a s (fib) is evaluated at 4.6 GeV where the bag parameters are provided as 
BX LL {m b ) = 0.87 ± 0.05, Bj_ R (m b ) = 1.75 ± 0.21, and B% R (m b ) = 1.16 ± 0.07 [22 1. 
(m s (2GeV) = lOOl^g MeV (l8| is evaluated as tti s (4.6 GeV) = MeV.) With these 

inputs, we obtain 



h s = (7.53 x 10 5 ) ■ (^) 2 + {gfb} 2 ~ k ■ g L sh g. 



L„R 
sb 



(34) 



where k = 5.05 ± 0.47. 

The value of h s should be as small as < 2.3 to satisfy the experimental constraint of ([TBI) . 
Therefore, the terms inside the squared bracket of f[34"j) must be as small as < 3 x 10~ 5 . This 
result can be rewritten as 



sb 



< 3.06 x 10" 6 . 



(35) 



Eq. fl3"5l) describes a complex hyperbolic surface which is flipped along the asymptotic complex 
lines satisfying 



te) 2 + {9sb) 2 -k-9s b 9sb = 



(36) 
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or equivalently, 



9sb = "fa , gf b = {l/a)g L sb , (37) 

where a = 4.84 for k = 5.05. On these asymptotic lines, 6l = Or where 6l,r is the phase 
of g^ R respectively. Consequently, the bound (135]) indicate that the generic values of \g^ R \ 
must be smaller than 10 -3 unless they are within (or close to) the asymptotic lines fl37j) . 
Since the AM S constraint parametrized by (1351) is highly stringent for \g^ R \ > 10~ 3 , the 
parameter space containing such values of couplings cannot avoid the fine tuning. 

Li Fl 

For the case that one of g sb ' is turned off, we can easily induce that the absolute value 
of the remaining nonzero coupling must be definitely smaller than 1.75 x 10~ 3 . Therefore, 
the required value of \g~l R \ or \g^ R \ for the explanation of the dimuon charge asymmetry in 



this case must be larger than 1, which is easily induced from analyzing the results in 



HQ. 



B. Ar s and from B s -> J/ip 

The enhancement of the like-sign dimuon charge asymmetry is constrained by the exper- 
imental measurement of the width difference AT S of the mass eigenstate B° s mesons, and 
the phase difference <p J J^^ between the B s mixing and the b — > sec decay. These are simul- 
taneously determined by measuring the indirect CP asymmetry of B s — y J/ifj <p decay. The 



recent result from the LHCb of lfb 1 integrated luminosity shows that 23] 



Ar s = 0.116 ± 0.018(stat.) ± 0.006(syst.) ps" 1 , (38) 
= _ .001 ± 0.101(stat.) ± 0.027(syst.) rad , (39) 

in which Ar s has about 1.2a deviation^ from (Ar s ) SM = (0.087 ± 0.021) ps" 1 and 0^ 
agrees well with the SM prediction (0^*) S m = Arg. ((V ts V t * b ) 2 / \V cs V* b ) 2 ) = -2/3f M = 
—0.036 ± 0.002 3]. Such new LHCb results dramatically reduce the room of new physics 
contribution in B s — B s mixing compared to those of the previous LHCb (337 pb _1 ), the 
CDF (5.2 fb- 1 ), and the DO (8.0 ft- 1 ). 

We first deal with the issue related with <pi , whose measurement at the LHCb lfb -1 
shows the most dramatic changes compared to the previous ones. The analytic expression of 



Note that the sign of Ar s is fixed to be positive in this result. 
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24j and |5|. Neglecting the SM strong phases in the B s — > J ftp <p 



4> J J^ ^ is well summarized in 
process, we obtain 0, Q 

sin<^ = sin(-2& + <j> a M ) + 2|r A | cos(-2/3 s + <j>' M ) sin^ A , (40) 

where (p s M = Arg.(Mi2/Mf 2 M ) is from the NP contribution in the dispersive part of B s — B s 
mixing, and the term with r\ is from the NP contribution in the b — > sec decay. Note tha t 
this result is obtained using the approximation that \r\\ <C 1 from the exact relation in 
In the figures to show the allowed parameter space, we use the exact relation. 



24). 



When there is no NP phase contribution in b — > sec process like our g TT scenario in Z' , we 

s 

Ad- 



have r\ = 0. Then, the NP effect in (f> J J^^ contributes only through <j> s M . Since we know that 



sin^j = h s sin 2o a j a/1 + h 2 s + 2h s cos 2o~ s = {h s sin2ag)y A ^ ~ /i s sin2cr s by combining 
the experimental bound of AM,, the value of h s sin2a s must be small enough to satisfy the 
measured result ( 1391) . In Fig. [3] (a), this result is represented in terms of our parameters 
2cr s and h s with the purple color surrounded by the thick red line. To satisfy both of the 
experimental bounds AM, f )T6|) and (j) J J^^ fl39|) . the allowed parameter space must satisfy 
h s < 0.3 1 except the small region around 2a s ~ it with 1.7 < h s < 2.2. 

Without NP contribution to Ti2, we can see that the enhancement in a s s£ is quite limited 
as — a s s i/ {a s s i) SM < 40 in Fig. |3](b), which is not satisfactory to explain the dimuon charge 
asymmetry within la. Considering the coefficient cot0^ M of h s in ( 1201) . we hence need at 
least (9(0.1) contribution by |rf 2 p /r^| to a s s£ . Applying our later result f l58l) in case of the 
g TT scenario, the value of coupling \g^ R g^ R \ should be much larger than 10~ 3 . Restricting 
|<7;^ fl | < 1 to avoid the nonperturbativity bound, we need \g^ R \ 3> 10 -3 which must be 
around the asymptotic lines (1371) to satisfy the AM, bound. Therefore, we can generically 
set 9l = Or in the g TT scenario, which induces 2d s = 2a s + nir for an integer n neglecting the 
contribution by the SM phases. We can conclude from this, at the region around 2a s ~ it 
with 1.7 < h s < 2.2, the enhancement in a s si = lm(Tl 2 /Ml 2 ) is ignorable in our g TT scenario. 
In result, we will proceed the analysis with the condition h s < 0.3 in this scenario. 

On the other hand, for the cases that we have NP phase contribution in b — > sec process 
like our g cc scenario, the NP effect in <p J s ^^ contributes also through r A ^ 0. The NP 
contribution in the B s — > J/ip amplitude is parametrized as 

J]((J/V0)a|Onp|^) = bxe** , (41) 
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(a) Experimental bounds in the B s sysmtem 



(b) Possible enhancement if Tf 2 = r^f M 



FIG. 3. The experimental bounds apply for the general NP scenarios without new phases in 
b — > sec, such as the g TT scenario, (a) The final allowed region of h s = \M^ /Mf^-\ and 2a s = 
Arg. (Mjf7 M l2 M ) from the AM *> (LHCb Tib" 1 , 0.371b- 1 , and DO Sftr 1 ) constraints is 

shown as the purple color surrounded by the thick red line. The yellow region (inside the dashed 
line boundary) : allowed by 90% AM S . The green region (inside the dot-dashed line boundary) 
: allowed by 90% (j) J s /i " p at the 8.0 fb" 1 DO. (0,^* = 0.15 ± 0.18(stat.) ± 0.06(syst.) [25j). The 
blue region (inside the line boundary) : allowed by 90% c^J^^ at the recent 1.0 fb _1 LHCb and 
the boundary at the 0.37 fb _1 at the LHCb in the last year is denoted as the dotted lines. The 
purple region surrounded by the thick red line denotes the allowed parameter space from all the 
commented constraints. The mainly remained region is h s < 0.3 <C 1, which provides a fine tuning 
choice in the parameter space explaining a s si . The other region of 2cr s ~ n with 1.7 < h s < 2.2 
is irrelevant in the enhancement of a s s » = Im(r| 2 /-^i2) i n the g TT scenario, (b) Without NP 
contribution to Ti2, we represent the possible enhancement of "(Jj/fe) with the numbers and 
the contours. In this case, the enhancement is quite limited such that — a s si /(a s se ) SM < 40 to be 
consistent with all the experimental bounds. 



where A is the polarization of final state vector particles. The longitudinal direction is A = 
and the two transverse directions are A = {+,—}• The angle <p\ is the new weak phase from 
the above NP contribution. The ratio of the amplitude \r\\ = b\/a\ is defined for the SM 
amplitude a\. 
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In the g cc scenario of Z' model, we obtain the following result according to 24| such that 



\r\=o\ 



1 M 2 z 2(gi + gX)(gL-k g 



isbJ 



where k = 1, k 
Melikhov-Stech 



8.8,9.8 and k 



9l 


Mf, 


V cb V* ■ 0.17 




1 


Mf 2( 




-9%) 


gl 


M\, 


V cb V* ■ 0.17 




l 


Mf 2( 


9cc + 9cc)(9sb- k - 


-9?b) 



9{Ml 



VrhV* ■ 0.17 



(42) 
(43) 
(44) 



= 0.11, 0.10 depending on the model of the form factors 
26j and Ball-Zwicky j^J , respectively! 7 ! The vector interaction of the charm 
quark pair is obtained from the factorization ( J /ip\crf^c\Q). 

Consequently, we obtain the following expression in the g cc scenario neglecting the SM 
prediction for \r\\ < 1. 

L sin2o"o 



sin 



a/1 + h\ + 2h s cos 2a, 
+ 



2 M*2(gg + g£)(gL b -k , ± gX) 



VrhV* ■ 0.17 



1 + h s cos 2a s 



(45) 



y/l + h% + 2h a COs2(X 



9\ M\, 

If we simply assume h s « which is conservatively safe from the AM S bound, we can 
neglect the first term of Eq. (j4"5"l) and the expression is simplified as 



sin 



(46) 



(1.0 x 10 3 ) \{g L cc + - fco i± sf 6 )| (sin^o,±) • 

To satisfy the recent LHCb result of lfb -1 with 90% C.L., we obtain the following simple 
condition on the couplings in this case 



1.7 x 10- 4 < \{g L cc + g*){g% - fc ,±^)| (sin^ 0j ±) < 1.7 X 10 



(47) 



which provides a strong constraint on the values of \g^ R gc2 R \ unless the Z' vector coupling 
to the charm quark pair is axial. For |g^J <C k + \gf b \, the most stringent bound is obtained 
from the A = + case and | siny9 + | « | sin^|. For the other case, the most stringent bound 
is obtained from the A = — case and | sin<^_| f» | svn9i\. Without considering the (almost) 
axial vector-like interaction of Z'cc, the constraint Q47D provides 



l^^sin^l <O(10- 5 ) 
|^ c L /sin^| < lO" 4 . 



(48) 



7 Actually, there are typically about 10 % theoretical uncertainties in the form factors. Such consideration 
in k+ as an example is shown in our Appendix [Bj 
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When 9l = Or or one of the couplings \g^ b \ and \gf b \ is dominant, the angle |sin(2cr s )| ~ 
| sm9i\ or | sui6(r|. In this case, hence, we can directly use the constraint (j4"Tj) to check the 
allowed parameter space for the dimuon charge asymmetry. 

In the mean time, we can also analyze more general case that the simple assumption 
h s ~ is not applied, while h s should still satisfy the AM S bound as Fig. [2j Then, the 
NP contribution in (j) J s ^ is small when there is a fine cancellation between the first and 
second terms in (j45p . From (1351) . we know that the off-diagonal couplings g^ b ,R must be 
around the asymptotic lines (13"T|) to satisfy AM S bound unless both of them are smaller 
than 10 -3 . Since the condition (13T1) demands 0l = 9r which makes the various constraints 
simpler, we can proceed our analysis according to the values of the off-diagonal couplings. 
For the clear readability of our paper, we leave the detail explanation in our Appendix O 
One thing to stress is that our numerical analysis in the g cc scenario will be proceeded with 
the conservative assumption h s ~ but our result can be generally applied even when a fine 
cancellation between the first and second terms in (1451) exists. 

Now, we move to the issue of Ar s . The analytic expression of Ar s /(Ar s ) SM is given as 

Ar s 2iri 2 |cos0 s 



(Ar s ) SM 2|rs 2 M |cos0fi 

i 



'1 + h s cos2a s )(l + h s cos2cr s ) + h s h s sin2cr s sin 2(7, 



a/1 + h 2 s + 2h s cos2ct s 

— tan (ji s sin 2o~ s ( 1 + h s cos 2d s ) — h s sin 2o s ( 1 + h s cos 2a s 



(49) 



With Eq. (|2Q|) . we can see that the enhancement of dimuon charge asymmetry is always 
possible without suffering from the constraint on Ar s /(Ar s ) SM . This is because the en- 
hancement of a q e is from Im(Fi2) and that of Ar g from Re(Fi2) along the direction of 
Re(M 12 ), as easily expected from the first relation in ( TH)]) . The consistent parameter space 
is shown with 2D plot as our Fig. HJ where the parameter space is free from the AM S bound. 
On the other hand, the other constraints from 28] such as B + — > K + t + t~, B s — > 



t+t~, B ->• X s t + t-, B X s7 , B ->- XJ + i~, and B K^*H + £- provide additional 
interesting limit in the allowed parameter space. (Among them, the strongest bound is 
given by B + — > K + t + t~.) The experimental bounds can be analytically expressed with 
h s cos2a s and h s sm2a s , in addition to the dimuon charge asymmetry value in (157|) . The 



bound is |r^ 2 NP /rf 2 SM | < 0.3 from 



Therefore, it is possible to check the consistency of 
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FIG. 4. This shows the consistency of explaining the observed a 1 L within la. The orange color 
with thick line boundary region denotes — aS s i/{o, s si ) SM > 580 and the dashed line 68 when h s = 
\Mf 2 p /Mf^\ <C 1 such as the g TT scenario. The experimental results with 90%C.L. AT S at the 



23] is 



LHCb lfb -1 is shown with light green color, while the result of B + — > K + t + t~ from 
shown in the light purple color. Actually this bound also covers other bounds such as B s — > t + t~ , 
B —7- X S T + T~ , B -> X s j, B -»• X s £+^", and S -> K^l + T. (A mong them, the strongest bound is 



given by B + — > K + t + t as seen in 



281].) Figure is simply depicted in terms of Re(r*2 /rf^ ) 



h s cos2a s and Im(r^ p /r*2 M ) = h s s'm2a s with the assumption that h s <C 1. We can easily see 



that the explanation of a s st and the experimental bound Ar s are orthogonal since they depend 
on the Im.(r^ 2 p ) and Re.(rj2 P ), respectively. The three experimental results are only marginally 
consistent. 



a s s£ , AT S , and B + — > K + t + t~ in terms of such parameters as Fig. |H Consequently, the 
three experimental results are only marginally consistent at the region allowing large NP 
contribution in a d st 
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FIG. 5. The colored region denotes the parameter space allowed by the 90% C.L. experimental 
bounds of 9uv\g^ R \- This is the case that 9l = Or = 6 = it/ 4 and g vv > 0. This parameter space 
is free from the Ar s bound by making rfl^/Tff 1 ^ almost imaginary, as well as the AM S bound. 
Since g^ T = g^u, the rough upper limit of the coupling is obtained g^ T \g^ R \ <3x 10~ 4 . Even for 
the other cases, the upper limit is below 10 -3 . 



svv 



L R T 

In the case that the non-zero g s ^ g^ T provides the sizable enhancement in a s sl) the coupling 
is constrained by its partner in the SU(2) doublet g„ u = g w . Following the analysis in 



9r 



29 



30J, we can obtain the limit of g V v9 s b from B — > K*vT>, B — > Kuu, and B — » X s vv. 



The detail way of calculating the Z' contribution in these processes are well summarized in 
our Appendix |D] 



As experimental upper bounds at 90% C.L. (1.65er), we obtain from 18J, |3lJ such that 



Bt(B -> JCvv) < 8 x 10~ 5 [18[ , 



Bt(B -> Kvv) < 1.3 x 1(T 5 
Br{B -> X„vv) < 6.4 x 10~ 4 (ai| • 



18| , 



(50) 
(51) 
(52) 



Combining all the limits, we obtain the limit of the couplings as Fig. [5j The allowed range 
of g vv \g^ R \ from the 90% C.L. experimental bounds is shown. We deal with the case that 
®l = 6r = ^ which is considered in the fine-tuned region of ( 157)) . This figure is an example 
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(a) Possible enhancement in a s sl (b) Limit of g^ T /g^ T 

FIG. 6. The limit of the couplings explaining the observed a% in the g TT scenario is shown. In 
(a), the region inside the blue line box is those remained after applying the b —> svv constraint, 
which is precisely shown in (b). In figure (b), we expressed the conservative exclusion region 
(grey region) based on the experimental bounds of 90% C.L from the b —> svv processes, which is 
\9tt9sI, R \ < 3 x 1CP 4 . This is the fine-tuned case that gf b = ag^ b and 9l = Or = tt/4 as an example. 
This parameter space is free from the Ar s bound by making rf^ p /rf| M almost imaginary, as well 
as the AM S bound. The ratio of ~ a ip/( a si) ^ s shown as a contour plot with the contours 68, 580, 
and 950. To explain the asymmetry within la, the value of |<jf^-| must be much larger than \g^ T \. 

9 = 7r /4 and g vv > 0. The rough upper limit of the coupling is obtained g vu \g^ b R \ < 3 x 10~ 4 . 
Even for the other cases, the upper limit is below 10~ 3 . 

In the g^ T scenario, this provides a strong direct upper bound of the couplings as shown 
in Fig. [6j To explain the asymmetry within la, the value of \g^ T \ must be much larger than 
that of |s£J. S 



D. sin 2/3 from B° J/i/jK s 



In this section, we deal with the additional experimental bound when the NP phases con- 
tribute to the b — > sec process, such as the g cc scenario. This is the indirect CP asymmetry 
8 In this case, the anomaly cancellation in the g TT scenario is threaten, unless we assume a scenario like 



the effective Z' model 
g? T coupling. 



32j | . This is because there is no way to cancel the SU(2) 2 U(1)' anomaly from the 
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sin 2/3 in the "golden plate" mode B — > J/ipKs- The SM prediction of sin 2/3 can be obtained 



from the fit of the unitarity triangle. According to [33j . we obtain sin(2/3) fit = 0.731 ±0.038, 
while the experimental measurements provide sin2/3 meas = 0.668 ± 0.028. In this case, the 
SM prediction is within la of the measured value. The detail analytic form of sin 2/3 is well 
described in {24] and jsj, which is similar to sin2/3 s in B s — > J/ip as (140 p . In the absence 
of the SM strong phase, 

sin 2/3 meas = sin(2/3) fit + 2|r| cos(2/3) fit sin <p . (53) 



As ( 140]) . this relation is obtained when |r| <C 1 and we use the exact relation in 
figures. 

In the g cc scenario, the analytic form of |r| is obtained as 

'sbj 



24J in our 



1 M 2 z 2(gg + g?M b + g^ 



, » - , - (5-2 x 10 2 ) x \(g^ c + g? c )(g^ b + g«)\ , (54) 

and the angle <p is simply obtained in the fine-tuned case (1371) such that <p = 9,9 + n. 
Therefore, the allowed range with 90% C.L. of the experimental result and the SM fit is 
obtained as 

-1.4 x 10- 4 < \{g L cc + flSfoJ + g*)\ sin <p < 1.4 x 10" 5 . (55) 

As the experimental bound by <pi , this provides the strong constraint on the NP param- 
eter space unless the coupling Z'cc is (almost) axial vector-like. This bound will be shown 
in Sec. |V] with other experimental constraints. 



On the other hand, the fitting value of sin(2/3) fit is enlarged if we drop the value of \V t 



Lib 



as an input since its inclusive and exclusive determination has a large difference. Instead, 
it is possible to use as inputs from the experiments, ex, AM s /AM rf , Br. (B — > tu). In this 
case, we obtain sin (2/3) fit = 0.867 ± 0.048 which induces more than 3a deviation from the 



observed central value [34| . By doing this, we can accommodate sizable NP contribution to 
sin 2/3 by g^ R gc£ R without sizable deviations in the B — > tv branching ratio and ex • In 
this case, the value of |r| from the NP contribution is allowed up to 20.0 + 6.5 = 26.5 % 
with the la predictions. In terms of the g cc scenario, the allowed range with 90% C.L. of 



9 In contrast to this interesting approach, it is fair to note that the Belle collaboration recently updated their 
result on Br(B~ — > t~v t ) which is consistent with the usual global fit to the Cabbibo-Kobayashi-Maskawa 
matrix elements 35]. 
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the experimental result and the SM fit in this case induces 



-2.8 x 1(T 4 < \(g L cc + g R )(g L sb + &)\wiip < -1.0 x 10" 4 . (56) 
The corresponding parameter region will be discussed in Sec. |V] 

IV. g TT SCENARIO FOR THE DIMUON CHARGE ASYMMETRY 

In this section, we explore the possible parameter space of the g TT scenario to explain 
the like-sign dimuon charge asymmetry, combined with the experimental bounds discussed 
in the previous section. In this scenario, the enhancement of T\ 2 is realized in the process 
of the t loop- induced Z' exchange. As seen in Fig. [5J we can simply assume h s <C 1 for the 
rough analysis. Then, the ratio of the flavor specific asymmetry 

« -(2.6 x lCr) h s sin 2a s + 1 + h s cos 2a s , 



where we put the central value of <f)^ M = 3.8 x 10 -3 . Due to the strong LHCb constraint on 
Ar s , the term —(2.6 x 10 2 ) h s sin 25" is dominant so that sin2cx s is far from 0. The value of 
h s in the g TT scenario is obtained from |5] such that 

h s « (6.7 x 10 3 ) x Abs. [{{g^f + (g R b f) " 0.5((^) a + (<£) 2 )} 

(58) 

{-3-3^ + 1.0((^ T ) 2 + (&)*)}] . 

To obtain the enhancement in a s si as large as our second reference point (1, —580) in 
fllOp . the rough lower limit of the couplings inside the Abs. symbol must be 3.3 x 10~ 4 = 
(1.8 x 10~ 2 ) 2 to for sin2cr s = 1. To obtain the enhancement as our first reference point 
(a d J(a d si ) SM ,a s J(a sl ) SM ) = (21,-68), the limit lowers to 3.9 x 10~ 5 = (6.2 x 10~ 3 ) 2 . Con- 
sequently, we roughly obtain the limit of the dominant new coupling to explain the dimuon 
charge asymmetry within la in the x 2_ fit 

\g^f\ > 1.8 x 10- 2 without (4) NP , 

(59) 

!U"<I IM > 6.2 x 10- 3 with 4/(4) SM = 21 . 

Since \g^ R \ < 1, the values of Max{|g^' fi |} cannot be smaller than ~ 6.2 x 10~ 3 so that 
the couplings g^ R lie on the asymptotic lines ( 157]) . having more than 1% fine tuning. The 
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FIG. 7. We changed the parameters in Fig. H]in terms of 2a s — \g^ b g^ T \ in the conservative case 
that g^ T = 0.1g^ T and gf b = ag^ b (h s « 0). The description on the colored region is same as that 
in Fig. [H This shows better understanding on the limits of the couplings in the g TT scenario. We 
see that the rough consistent region of \gf b g^ T \ is about 10~ 2 with a^ e /(a^ e ) SM = 21. 



allowed parameter space is shown in Fig. [7] where we used g^ T = 0-lg^ T to maximally satisfy 
the constraint by b — > svv as explained in Sec. IIII CI We see that the rough consistent region 
of \gf b g^ T \ is about 10~ 2 with G^ £ /(af £ ) SM = 21. We also show the allowed parameter space 



of Ar, and B + ->- K+t+t' from 



redrawn from the allowed region in Fig. HJ 



Consequently, the g TT scenario where the Z' coupling to the r pair enhances the a s sl 
requires the existence of the coupling \g^ b R g^ R \ larger than about 10~ 2 to explain the dimuon 
charge asymmetry. Therefore, this parameter space cannot avoid the fine tuning from the 
AM S constraint. In addition, due to the constraint from the b — > svv experiments, the 
coupling \g^ T \ must be as small as 3 x 10~ 4 . This result demands a non-trivial approach 
in establishing an anomaly free model as mentioned at the end of Sec. IIII CI The allowed 
parameter space explaining the observed a s s£ within la is marginally consistent with the 



experimental bounds of Ar s at the LHCb lfb 1 and the result of B + — > K + t + t from 28 ]. 
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FIG. 8. We represent the allowed parameter space explaining the dimuon charge asymmetry by 
fixing g^ c = and g^ c > on the fine tuning region gf b = ag^ b . (Therefore, 9l = Or = 6.) 
The parameter region g^ c \g^ b \ > 10 -2 is not considered to avoid the rough constraint from B° — > 



D + D~ [38]. Even though we allow the NP contribution which is about half of the SM tree level 
prediction, it is roughly g^ c \gf b \/g\ < 0.5\V c i,V* s \ ~ 0.02. The numbers in the contours denote the 
ratio — flf£/(a^) SM as previous figures. The light pink region denotes the 90% bound from the 
B° -)• J/ipKs considering the usual fit sin(2/3) fit = 0.731 ± 0.038 and the area surrounded by the 



341 ]. The blue region is 



red dashed line is for the special fit sin(2/3) fit = 0.867 ± 0.048 in 
of (j) J J^^ at the recent 1.0 fb" 1 LHCb. The minimum value of {g^g^l to explain the DO dimuon 
charge asymmetry within la is about 8 x 10~ 3 and 10~ 3 , without the NP contribution to a?p and 
with the maximal contribution of «^/(o^) SM = 21, respectively. These bounds do not satisfy the 
experimental constraints, which is expected in our simple analysis. 



V. g cc SCENARIO 



In this section, we explore the possible parameter space of the g cc scenario to explain the 
like-sign dimuon charge asymmetry, combined with the experimental bounds discussed in 
the Sec. IIHI The enhancement of T* 2 from the interference of the SM process and the Z' 
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induced tree level FCNC in b — > sec is calculated from [361 |37J such that 

ml „, 



3n(2m Ba ) g{ M|, 



x 



(60) 



444 {(1 " Xc)(O ll ) + (1 + 2x c )(O fifl )} 
+ igXgg {(1 - a: c )(0 Li? ) + (1 + 2i c )(O m >} 
+12a^£<<^> + 12x c 44(0 LK >] , 

where x c = m 2 /m 2 and .fTi = 3.11 is calculated from the Wilson coefficient of the corre- 
sponding operators as in jsfl This result is different from that j^] especially adding the 
suppression factor x c at the coefficient of the contribution by 44- 
The value of h s in the g cc scenario is obtained such that 

h s « 173.7 x | (I.15^ b - 1.764) 4 + (-1-034 + 0.644) 4| . (61) 
When /i s <l, the ratio of the flavor specific asymmetry is 

<i/Kif M « -4-6 x 10 4 x | (1.154 - 1.764) 4 + (-1-034 + 0.644) g L cc \ sin2a s , 

(62) 

in the region that the contribution from Re(r^ p /r^| M ) is suppressed to avoid the Ar s 



bound. Then, we directly obtain the rough lower limit of the couplings from f[62|) . To obtain 
the dimuon charge asymmetry within ler without any NP contribution in a d e , we have 

| (1.154 - 1-764) 91 + (-1-034 + 0.644) g L cc \ • | sm9 L{R) \ > 1.3 x 10~ 2 , (63) 
while for (a d si /(a d e ) SM ,a s J(a sl ) SM ) = (21, -68), we have 

| (1.154 - 1-764) 4 + (-1-034 + 0.644) 4| • | sm6 L(R) \ > 1.4 x 1(T 3 . (64) 

Consequently, we roughly obtain the limit of the dominant new coupling to explain the 
dimuon charge asymmetry within la in the x 2 -fit 

14*4* sin e L[R) | > 10- 2 without (4) NP , 

(65) 

\g^ R g^ R sm6 L{R) \ > 1(T 3 with a d J (a d s£ ) SM = 21 . 

In this case, the magnitude of 14*1 nas ^ ess ^ ne tuning from the AM, constraint compared 
to that in the g TT scenario. Unless the interaction Z'cc is (almost) axial vector-like, the 
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FIG. 9. We show what extent the interaction Z'cc should be axial vector-like in this figure. For 
various values of the difference 5 C = (g^ c + g^ c )/g^ c , our parameter space explaining the dimoun 
charge asymmetry is shown for 5 C > 0, by fixing g^ b ~ (l/a)gh with 6l = Or = 3ir/2, which 
is different from the case in Fig. [8j The numbers in the contours denote the value of the ratio 
— a^/(a^) SM . The contour line with pink color is the boundary — a s si /(a s si ) SM = 68 which demand 
a te/( a te) SM = 21 to explain the dimoun charge asymmetry within la in the fit of Fig. [TJ The 
meshed blue region denotes the 90% allowed region by the recent LHCb lib -1 result of ^ 
and the area surrounded by the red dash line is that explaining the sin 2/3 with the special fit 

. Of course, all the blue region of <\/ s ^ is also allowed 



sin(2/3) fit = 0.867 ± 0.048 proposed in 



by that using the usual fit sin(2/3) fit = 0.731 ±0.038. In this case, the coupling \g^ c gf b \ > 7.5 x 10~ 3 
and 10 -3 are required for the explanation of the asymmetry within la without the NP contribution 
to orn and with the maximal contribution of ^»/(co) = 21, respectively. Our parameter space 
with — a^/( a s^) SM > 580 and a d sl /{a d s ^f M = 1 can simultaneously explain the recent LHCb lib -1 
result and the sin 2/3 for 5 C < 2.5 x 10~ 2 . 

result ( )65|) shows a direct contradiction with the constraint by <f> J s we obtained in ( l4~8j) . as 
well as the sin 2/3 measurements in Sec. IIIIDI This result is shown in Fig. [8] in the simple 

case g%. = and gf b = ag L sh . 

When \gc C +gcc\ "C 1, the constraint by the <pi^^ at the 1.0 fb _1 LHCb is loosen. We show 
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what extent the interaction Z'cc should be axial vector-like in our Fig. [9j simultaneously 



explaining the interesting parameter region of sin 2/3 using sin(2/3) fit = 0.867 ± 0.048 34]. 
As a result, for the explanation of the asymmetry within la without the NP contribution 
to a d sl , we can find the consistent parameter space {g^g^l > 7.5 x 10~ 3 and 5 C < 2.5 x 10 -2 
when g R b = ag^ b with 6 L = 6 R = 3%/2. 

Finally, we discuss the possibility of the model construction providing the (almost) axial 
vector-like interaction Z'cc. In this case, we need the sizable coupling \g^ c \ ~ \g R \ at least 
0(1O~ 2 ). Since the left-handed charm quark constitutes a SU(2) doublet with the left- 
handed strange quark, the coupling g^ s is also of C(10 -2 ). Then, we obtain the nonzero 
off-diagonal couplings g^ c and g% s , unless the size of the couplings gu U i.9dd) are same as 
9cc(9ss)- ft mus t be noted that the size of such off-diagonal couplings are constrained by 
the bounds such as K or D meson mixings. For example, the D meson mixing provides the 



strong upper limit of the g^ c coupling as small as 2 x 10~ 4 [3_9|, |40j. Considering the CKM 
relation of g^ c ~ 0.23g^ R , the coupling \g<£ R \ must be smaller than 10~ 3 , which in turn 
demands \g^ R gci R \ < 10~ 3 . This is out of the la region for the dimuon charge asymmetry 
even with arbitrary contribution in a d e , as seen in Fig. [8]and[9l 

On the other hand, we need to consider the bounds of the 7r production processes from 
the B s meson when the size of the couplings g£ u = g\ d are almost same as g^ c = g^ s . For 



example, the upper bound of Br(B s — > tt + tt ) is as strong as 1.2 x 10 6 18|, while the value 
of I 

9uudd 9^1 is as large as \g R gf b \ > 7.5 x 10 3 . Therefore, the scenario with the (almost) 
axial vector-like interaction Z'cc is not plausible. 

Consequently, the g cc scenario where the Z' coupling to the c-quark pair enhances the a s s£ 
requires the existence of the coupling \g^ R g^i, R \ larger than 0(1O -3 ) to explain the dimuon 
charge asymmetry. This parameter space has smaller fine tuning from the AM S compared 
to the g TT scenario, due to the interference with the SM process in the contribution to T* 2 . 
However, the recent LHCb lfb -1 constraint on <pl^^ and Ar s , as well as the constraint from 
B — > J/ipKs, is quite strong to demand the (almost) axial vector-like interaction of Z'cc. 
On the other hand, the existence of such interaction makes the model construction very 
hard due to the experimental bounds such as K or D meson mixing and the 7r production 
from the B s decays. 
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VI. CONCLUSIONS 



The like-sign dimuon charge asymmetry has been observed at the DO which is deviated 
more than 3a from the SM prediction. In the recent result in 2011, it was possible to 
separately detect the flavor specific asymmetry from the B s and mixing by imposing 
the impact parameter cut reducing the background. In this paper, we showed that the 
enhancement of flavor specific asymmetry a s s£ is highly constrained by the recent LHCb 
result with lfb -1 integrated luminosity. We presented the constraints on the Z' couplings 
9bs, 9tt or g cc and the possible enhancement of a s st The actual upper bound of the couplings 
are expressed when M z > ~ M z . By simple scaling of the ratio M z >/M z , our result can be 
applied to the other mass of Z' as well. 

For the flavor specific asymmetry a s s£ , there are three kinds of criteria. By allowing 
sizable new physics contribution in Bd system (a d e is 21 times larger than the SM prediction 
from NP), \a s s g/a s se SM \ > 68 is needed to be within la region. If there is no new physics 
in B d system and the deviation of A b s£ is the consequence of B s system alone, we need 
\a s se / a s se SM \ > 580 to be within la region. The central value requires \a s s g/a s se SM \ > 950. 

The B s system is highly constrained by the recent LHCb data. In the absence of the 
modification in the decay, T[ 2 , the recent LHCb measurement of <\> J J^^ strongly constrains 
the phase of the mixing, M( 2 . As a result, the maximum enhancement of a s si from M[ 2 alone 
is at most 40 times the SM prediction which is not enough to be within la even if we allow 
arbitrary NP contribution to a d sl . 

The b — > ccs coupling can provide an extra contribution in (f> J s ^^ from the decay. If 
the coupling is small enough, the main effect would be to modify the relation between the 
phase of M{ 2 and (j} J s ^^ and the constraint on (j} J s ^^ can be slightly relaxed. In this case the 
main impact of b — > ccs is to avoid the constraints on the phase of M( 2 from the LHCb 
measurement of 0^^. By allowing the b — > ccs coupling, a s si can be as large as 50 times the 
SM prediction, which is still smaller than 68 for the la explanation of the asymmetry with 
arbitrary a d e . 

The T s l2 is constrained by AT S measurement which is basically Re(r^ 2 ) when M{ 2 is 
almost real. The enhancement of a s s£ is mainly from Im(rf 2 ) which can affect the other 
observables like B + — > K + tt in the g TT scenario and B s — > J/ip in the g cc scenario. 

The g TT scenario where the Z' coupling to the r pair enhances the a s st requires the 
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existence of the coupling \g^ R g^ R \ larger than about 1CT 2 to explain the dimuon charge 
asymmetry. Therefore, this parameter space cannot avoid the fine tuning from the AM S 
constraint \g a j,\ < 10 -3 . In addition, due to the constraint from the b — > svv experiments, 
\g^ T \ must be as small as 3 x 10 -4 . The allowed parameter space explaining the observed 
a s sl within la (68 times larger than the SM prediction) is marginally consistent with the 
experimental bounds of Ar s at the LHCb lfb -1 . 

The g cc scenario where the Z' coupling to the c-quark pair enhances the a s sl requires the 
existence of the coupling \g^ R g^ R \ larger than about 1CT 3 to explain the dimuon charge 
asymmetry. This parameter space has smaller fine tuning from the AM S compared to 
the g TT scenario, due to the interference with the SM process in the contribution to T* 2 . 
However, the recent LHCb lfb -1 constraint on ^ and AT S , as well as the constraint from 
B° — > J/ipKs, are quite strong. So the interaction Z'cc must be (almost) axial vector-like. 
On the other hand, the existence of g^ c from the axial vector constraints makes the model 
construction not plausible due to the experimental bounds such as K or D meson mixing 
and the tc production from the B s decays since g^ c = g^ s . 

Consequently, it is impossible to explain the la range of like-sign dimuon charge asymme- 
try using Z' contribution in B s system without the enhancement in a d e . Even with arbitrary 
a d e , the g cc scenario demands unrealistic model construction. Therefore, we need to consider 
the sizable NP contribution in a d £ , while making the a s sl as small as possible. To explain 
the asymmetry within la by minimizing the NP contribution to B s system, we need the 
a d sl which is only about 21 times the SM prediction at most, as shown in Fig. [TJ So the 
required off-diagonal coupling \g^,\ to enhance the a d e is smaller than the \g s b\- On the other 
hand, the CKM suppression strengthens the experimental bounds except Ar^ which has 
been poorly measured so far. The experimental bounds to be analyzed contain B° — > t + t~ , 
B — > 7it + t~, and B° s — > K t + t~ when we consider enhancement in Tf 2 through nonzero 
g TT coupling. For the case with nonzero g cc coupling, we need to consider B — > J/ipn and 
B s — > KqJ/iP, etc. Therefore, more careful analysis is required in the future to investigate 
this approach. 
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Appendix A: Calculation of the contribution in to b — > 57 

Considering the non-zero value of g b £ R or <?f s ' fl , we can obtain the upper limit of \gf' Ll 



) Hi, 



from the b — > 37 penguin constraint, shown in Fig. [TUJ 




FIG. 10. b — > S7 penguin contribution 

The inclusive decay T(B — > X s 7) is given approximately by T(b — > X^ arton/ y). The 
nonperturbative correction to this approximation is smaller than the NNLO perturbative 
QCD corrections to T(b — > XJ arton 7). The theoretical prediction for the partonic T(b — > 
Xparton^ j g usua Uy normalized by the semileptonic decay rate to get rid of the uncertainties 
related to the CKM matrix elements and the fifth power of the 6-quark mass. Therefore, 



the SM NNLO result for a photon-energy cut of > 1.6 GeV is obtained as [42| 



Br(5 X s7 ) SM = (3.15 ± 0.23) x 10" 4 , (Al) 



while the experimental result with for the same energy cut is measured as 43] 



Br(5 ->• X s7 ) cxp = (3.55 ± 0.24 ± 0.09) x 10~ 4 . (A2) 



The NP contribution in the total branching ratio is below 30% seeing the result of (I A 1 1) 
and (lA2p . Therefore, a naive strongest constraint of \g^ R g^i, R \ is < 10~ 2 as shown in Fig. 
[TT1 However, larger values of the couplings can still satisfy the b — > 57 constraint once the 
coupling ratio gf h jg^ h is about 1.1 or 1.27. The SM NNLO contribution shows a negligible 
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FIG. 11. The limit of the couplings from the experimental bounds of 90% C.L (Blue) and 95% 
C.L. (Cyan, Dashed boundary line) of B — > X s ~f for fine-tuned cases (a) 6l = Or = vr/4 and (b) 
6 L = 6 R = 3tt/4. 



dependence on the value of \x b . When the LO NP contribution enhances the SM value 
by 20%, the /x& dependence in the total branching ratio induces about 3% uncertainty for 
fib = 2.5 — 5 GeVin the numerical analysis in [441 ]. 



Following Eq. (5.3) of 44j, the branching ratio with the NP contribution is obtained 



Bi(B X s7 ) = (2.47 x 10~ 3 ) 

x (\C 7 ^ b )\ 2 + \C' 77 (fi b )\ 2 + N(E 7 )) 



(A3) 



where N(E 7 ) = (3.6 ± 0.6) x 10 3 is a nonperturbative contribution. Considering the LO 
NP contributions 



C 77 M = C 7 s 7 M (/i 6 ) + AC 77 (/i 6 ) 



(A4) 



where the central value of the SM contribution is calculated at the NNLO level for fif, 
2.5 GeV such that 



Cf-fM = -0.3525 . 



(A5) 



30 



The NP contribution is obtained as following. 

AC 77 (/i fe ) 



1 Mf 1 



+ ) 9ss9sb - ^ S LL9ss(9s b Y 



2 \ m 

+ ( -« 7 - 2k 8 — - 2/4^0^) 



~/^7 — 2/xg — 2/x 



Li? 



sb) 



where k's are listed in Table 1 of 



44j |. The prime coefficients are obtained as 



(A6) 



O^) = -0-3523 



(A7) 



2m, 



(-«r + 3k 8 ) ^« - 2<^f s (^)* 
-«r + 3/%) - 24^ ^fo*)* 



9m h 



+ 



2m s 
9m>, 



2 m. 



3 m 6 



m 6 



In the simple case that g£ R = and our values of ks are not much different from those with 



+ ( k 7 - 2—k 8 - 2k" lr ) gk(g*y + ( -k? - 2k 8 )( — ) gf s g^ - 2 K s LR gi(g 



m. 



isb) 



(A8) 



the matching scale at around 200 GeV. We obtain the following result, 



C 77 = -0.3523 - 9.11 x (gt - 1.27g*){g, 



sb) J 



C' 7l = -0.3523^ + 6.83 x - 1.1^)(4)* 



m b 



(A9) 
(A10) 



Plugging (1A9j) into (1A3|) with the consideration of the experimental limit (1A2|) . we can 
obtain the limit of (g bb — 1.27 g^)\g^ h \ and (g bb — l.lgH)\g^ b \ according to a fixed value of the 
Or) such that 



6.42(^ - 1.27^)|^ 6 | cos0 L + 82.99(<j& - 1.27^) 2 |^ b | 2 

- 0.11(flS - l.lflg)|fl2| cos0 fl + 46.65G/& - Llfl2) 2 |fl5| a = 0.016 ± 0.01 



(All) 



within la up to C(10~ 4 ), calculated with m s = 100 MeV and m b = 4.2 GeV. 
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Appendix B: The effect of the theoretical uncertainties in the form factors 



Considering the square root error propagation, this uncertainty changes the quantity k± 
to 1/3 or twice of the original calculation. (For the Ball-Zwicky model, k + = (8.02Ai + 
3.35Z) / (8.02Ai — 3.35Z) for the form factors A\, Z in |24|. Considering the 10 % uncertain- 
ties, k + = (8.02-0.42(l±0.1) + 3.35-0.82(l±0.1))/(8.02-0.42(l±0.1)-3.35-0.82(l±0.1)). 
Simple calculation with A\ = 0.42(1 — 0.9) and Z = 0.82(1 + 0.1) can make k + very large as 
518. However, the error propagation without considering the covariance can make it smaller. 

8.02 ■ 0.42 + 3.35 ■ 0.82 ± v/(8.02) 2 (0.042) 2 + (3.35) 2 (0.082) 2 



k, 



8.02 • 0.42 - 3.35 ■ 0.82 ± v/(8.02) 2 (0.042) 2 + (3.35) 2 (0.082) 2 
8.02-0.42 + 3.35-0.82 



8.02 ■ 0.42 - 3.35 ■ 0.82 



1 ± v / (8.02) 2 (0.042) 2 + (3.35) 2 (0.082) 2 




(Bl) 



8.02 • 0.42 + 3.35 • 0.82) 2 (8.02 • 0.42 - 3.35 • 0.82) 2 



= 9.83(1 ±0.70) , 

which makes k + = 2.9 — 16.7.) The ratio of each polarized amplitude can be obtained in the 
CDF data where the transverse amplitude A\\ t ± = (A+ + A-)/y/2 41]. 



Appendix C: The analysis on <j> s when the assumption h s « is not applied 

As commented in the article, we analyze more general case that the simple assumption 
/i s ~ is not applied, while h s should still satisfy the AM S bound as Fig. |2j Then, the NP 
contribution in § J J^^ is small when there is a fine cancellation between the first and second 
terms in ( I4"5j) . From (I3"5"|) . we know that the off-diagonal couplings g R £ R must be around 
the asymptotic lines (137)) to satisfy AM S bound unless both of them are smaller than 10~ 3 . 
Since the condition ( )37|) demands 6l = 6r which makes the various constraints simpler, we 
can proceed our analysis according to the values of the off-diagonal couplings. Therefore, 
we classify the cases as following. 

i) At least one of \g^ R \ > 10~ 3 but g^ R is small enough to ignore rf 2 p /r^. |^ 

ii) At least one of \g^ R \ > 10 -3 and g^ R is large. (rf 2 p /r^ is sizable.) 
10 The coupling \g^ R gcc R \ can be either small or are in a special relation making Ti2 small. 



32 





0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5 

,„ NP,., SM .., NP,., SM I 

(a) Case i) (b) Case iii) 

FIG. 12. We represent the allowed parameter space explaining the dimuon charge asymmetry 
within la for the case i) (a) and iii) (b) in the contents, (a) The purple region surrounded by the 
thick red line is the case without the new phase in b — > sec as explained in Fig. When h s is 
sizable to cancel the contribution from the b — > sec, larger region is allowed by (p J s l ^' t> but the NP 
parameter space is still constrained by sin 2/3 as explained in Sec. IIIIDI The blue line represents 
the combined bound of AM S , <pi^^, and sin 2/3 for the case i). The contours denote the 50, 68, 
580, and 950 of the ratio — a^/(a^) SM . Even with arbitrary contribution in a%, we see that the 
combined bound do not allow the enhancement — a s s g/(a 8 s g) SM = 68 to explain the dimuon charge 
asymmetry within la. (b) The green region surrounded by the dotted line is the 90% allowed 
region of Ar s from the LHCb result (|38p . We can see that the NP contribution in mixing is 
necessary to explain the dimuon charge asymmetry within la. 

iii) \g^ R \ < 10~ 3 and g^ R is small enough to ignore rf 2 p /r^. 

iv) \g^ R \ < 1CT 3 but g^ R is large enough to make rf 2 p /r^ sizable . 

For the case i), the off-diagonal couplings g^ R must satisfy the fine tuning condition (I57|) . 
which demands 8l = Or. Then, it is possible to directly apply the constraint of sin 2/3 from 
B° — > J/ipKs as explained in Sec. IIII Dl Considering this effect, it is possible to obtain the 
limit of the parameters 2a s and h s as shown in Fig. [12] (a). In the figure, the combined 
bound from AM S , <pi^^, and sin 2/3 is represented as the region surrounded by the blue line. 
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In result, it is impossible to obtain the enough enhancement of a s sl to explain the dimuon 
charge asymmetry within la in case i). 

For the cases ii) and iv), the value of h s takes the dominant role in enhancing the a s s£ . 
(h s must be at least 0(1).) Therefore, the order of couplings g^g^ m ust be larger 
than 10~ 3 as our numerical result in Sec. |V] For the case ii), the constraint from sin 2/5 
in Sec. IIII Dl excludes the parameter space of \g^ R 9cc R \ > C(10" 3 ) just as our result in 
Sec. [V] In addition, we do not need to consider the case iv) since \g^ R \ > 1 to satisfy 
\9% R 9% R \ > O(10-*). 

For the case hi), we cannot directly apply the constraint from sin 2(3 since the condition 
L m Or is not necessary. We now apply — a^/( a ^) SM > 68 to Eq. (fl5|) so that the value 
of -sin0^ = - Arg. (M 1 s 2 /Mf 2 SM ) must be of order Q8(AM s /AMf M ) sin 0f M > 0.2 from 
Eq. (HSJ since r^/r^ is negligible. If we demand a fine cancellation in (145 p . the value 
of sin cj) s M ~ h s siu2a s must be similar to (1.0 x 103)\(g^ c -\- g R c )(g^ b — kxg^) sunpx\, which 
i s ^ Kficc + 9 R ) s ^ n fx\- Therefore, one of the couplings \g^ R \ must be at least of order 
(9(10 _1 ). Then, we see that the value of h s is as small as O(10~~ 2 ) from Eq. ( 16T|) . Now, 
simply assuming h s w 0, we can directly compare the experimental bounds of AM S and 
AT S analytically. As seen in Fig. [12] (b), the NP contribution in B^ mixing is necessary 
to explain the dimuon charge asymmetry within la. This result is same as ours in the 
numerical analysis with the assumption h s « in Sec. |V] 



Appendix D: Calculation of the contribution to b — > svv 

All the observables depend on the complex Wilson coefficient Cl and Cr such that 

Br(fi ->- K*vv) = 6.8 x KT 6 (1 + 1.317?) e 2 , (Dl) 
Bi(B ->- Kvv) = 4.5 x 10~ 6 (1 - 2rj) e 2 , (D2) 
Br(B X s vv) = 2.7 x 10" 5 (1 + 0.097?) e 2 , (D3) 



where e and rj are 



V\c L \ 


2 + l< 


C R \ 2 


\(C L ) SM 





(D4) 



_ -Re(C L C R ) 

v ~ \c L \* + \c R \* ■ (U5) 
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The Wilson coefficients are 



1 1 M 



2 



C L = (C L ) SM + (C L r - (^) SM - t^±—±^g w A , (D6) 

1 2~^ V tb V ts 9l M Z' 



1 1 1 M Z „ „R 



Cr - _ o a t/ t/* ~Z2T7T 9™9sb ) ( D7 ) 

where (Cl) sm = —6.83 ± 0.06 It is clearly seen that the SM prediction is obtained when 
7] = and e = 1. 
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